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Introduction 


•  In  many  engineering  applications  involving  fluid-thermal  systems,  detailed  quantitative  infomation 
on  the  flow,  temperature  and  species  concentration  is  needed  for  system  optimization 
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Introduction 


•  Numerical  simulation  can  obtain  the  desired  information  and  thus  optimize  the  system 

However,  this  approach  requires  well-defined  boundary  and  operating  conditions  which  may  not 
be  completely  known  due  to  limited  access  for  experimental  measurements 
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Introduction 


•  The  objective  of  our  research  is  to  develop  a  Dynamic  Data  Driven  Applications  System 
approach  that  synergizes  experiment  and  simulation  to  determine  the  boundary  and  operating 
conditions,  thereby  achieving  a  full  simulation  capability 
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Problem  Definition 


Jet  in  Crossflow 

•  Heated  wall  jet  in  crossflow 

The  objective  is  to  determine  the  Jet  inflow  conditions  {Uj,  Tj)  using  a  Dynamic  Data  Driven 
Applications  Systems  method  that  synergizes  experiment  and  simulation 
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Problem  Definition 


Jet  in  Crossflow 


•  Experiment 

Rutgers  Low  Speed  Wind  Tunnel 
Non-intrusive  laser  diode  measurement 
Measure  absorbance  \/s  time  at  fixed  (x,y) 

Static  temperature  T  \/s  time  from  absorbance 
Limited  region  for  absorbance  measurement 
Each  y)  measurement  reguires  1  hr  Experimental  configuration 
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Problem  Definition 


Jet  in  Crossflow 


•  Laser  diode  absorbance 
Instantaneous  absorbance 

N  (lo-  I(x,y,t)) 

y)  = - - - 

J-o 

where  lo  is  incident  intensity  at  (x,y,zi)  and  I(x,y,t) 
is  transmitted  intensity  at  {x,y,z2) 

Absorbance  per  cm  of  the  ^i?2(6)  line  of  the  oxygen 
transition  biY.^ u'  =  0  ^  =  0  at  761.139  nm 

is 

—  =  0.083T“^  -  2.26  •  10“^ 

dz 

where  T{x,y,z,t)  is  the  static  temperature  in  K 


Laser  diode  arrangement 


Typical  T  vs  time 
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Problem  Definition 


Jet  in  Crossflow 


•  Simulation 

Laminar  Navier-Stokes  equations 
Incompressible,  ideal  gas 
Unsteady,  time-dependent 
Sutherland  viscosity  law 
Fluent© 

Parallel  (8  processors) 


E  D 
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Problem  Definition 


Jet  in  Crossflow 


•  Flow  Structure 


t  =  80  ms 


t  =  120  ms 
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Problem  Definition 


Jet  in  Crossflow 


•  Assumptions 


Large  set  Ss  of  discrete  data  locations 
defined  (<  no.  of  grid  cells  in  simulation)  e 


For  each  experiment,  time  series  data 
obtained  for  small  subset  =  1,2,... 
of  locations 

For  each  simulation,  time  series  data 

obtained  for  entire  set  Ss  for  each  Uj  and 

T 

-^3 

•  The  quantity  for  comparison  between 
experiment  and  simulation  is  the  mean 
temperature  Tm(x,y) 


Boundary  of  computational  domain 
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•  Problem 


Develop  and  apply  a  DDDAS  Methodology 
for  determining  Uj  and  Tj 
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Response  Surface  Models 


•  Energy  equation  decouples  from  the  mass  and  momentum  equations 

•  Instantaneous  temperature  behaves  as  passive  scalar  and  thus  must 
scale  as 

Flow  Conditions 


T(x,  y,  t)  -  Too  =  {Tj  -  Too)f(x,  y,  t]  Uj,  Uoo) 

•  Response  Surface  Model 


Tmi^X,  y)  Too 


l3o{x,y)  +  I3i(x,y) 


+  P2ix,y) 


•  The  coefficients  I3i(x,y)  are  obtained  from  simulations  performed  for 
a  fixed  value  T-  —  Too  (selected  from  the  range  indicated  in  Table) 
and  a  set  of  Uj 


Parameter  Value 


Uoo  (m/s) 
Too  (K) 
Poo  (kPa) 
Uj  (m/s) 
T,  (K) 

Pj  (kPa) 


4.0 
290. 
101.8 
4.0  to  8.0 
350  to  450 
101.8 
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Dynamic  Data  Driven  Applications  System  Methodology 


E  D 


1.  Select  monitor  locations  Ss  for 
simulations 

2.  Generate  Response  Surface  Models 
based  on  simulations  for  fixed  ATj 

3.  Select  monitor  locations  for 
experiments 

4.  Estimate  experimental  values  for 
Tj—Too  and  Uj  using  Response 
Surface  Models  and  experimental 
data  at  monitor  locations 

5.  Repeat  at  Step  No.  2  if  estimated 
Tj  —  Too  is  significantly  different  than 
used  to  generate  Response  Surface 
Models;  otherwise,  determine  new 
measurement  locations  Se~^'^ 

6.  Repeat  until  converged 


No. 

X 

y 

No. 

X 

y 

No. 

X 

y 

1 

1.2 

2.0 

7 

1.2 

3.0 

13 

1.2 

4.0 

2 

3.2 

2.0 

8 

3.2 

3.0 

14 

3.2 

4.0 

3 

5.2 

2.0 

9 

5.2 

3.0 

15 

5.2 

4.0 

4 

7.2 

2.0 

10 

7.2 

3.0 

16 

7.2 

4.0 

5 

9.2 

2.0 

11 

9.2 

3.0 

17 

9.2 

4.0 

6 

11.2 

2.0 

12 

11.2 

3.0 

18 

11.2 

4.0 

Distances  in  cm  from  jet  center 
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Dynamic  Data  Driven  Applications  System  Methodology 


•  Estimating  experimental  value  of  T-  —  Too  and  Uj 


—  Calculate  square  error  between  the  experimental  mean  temperature  and  the  Response  Surface 
Model  for  each  possible  subset  of  I  locations  within  Sg  as  computed  as 


E  =  J2{  ^Tme  -  at. 


j 


Po(x,y)  +  Pi(x,y)  I  +  P2(x,y)  ' 


a 


oo 


a 


oo 


where  ATj  =  Tj  —  Too,  ^Tme  =  Tme  —  Too,  and  the  sum  is  over  I  locations  within  Sg  (the 
minimum  number  for  I  is  2) 

Example:  Assume  contains  six  locations  and  let  I  =  2.  For  each  possible  set  of  two  locations  from  S^,  the 
values  of  AT^  and  Uj  that  minimize  E  are  determined.  This  yields  fifteen  triplets  E). 


—  For  a  given  value  of  the  predicted  values  of  ATj  and  Uj,  denoted  by  ATj  and  f/j,  are  taken 
to  be  the  triplet  with  the  minimum  E  (/.e.,  the  values  of  ATj  and  Uj  with  the  smallest  square 
error). 


—  The  procedure  is  repeated  for  all  values  of  I  from  /  =  2  to  n  =  size  S^. 


—  The  estimate  for  the  experimental  value  of  Tj  —  Too  is  the  average  of  these  values 
Tj  —  Too  =  (n— 1)“^  ATj  and  similarly  for  Uj. 
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Results 


•  Application  of  DDDAS  Methodology 

No.  Step  Description 

1  1  A  total  of  eighteen  monitor  locations  were  selected 

2  2  Response  Surface  Models  were  generated  at  all  monitor  locations  using  ATj  =  66  K 

3  3  Six  locations  (Nos.  3,  9,  10,  14,  15  and  16)  were  selected  for  experiment 

4  4  Using  the  experimental  mean  temperature  measurements  at  the  six  locations,  the 

estimated  values  ATj  =  110  zb  16  K  and  Uj  =  7.3  zb  1  m/s  obtained  using  the  RSMs 

5  5  A  new  set  of  locations  for  experiments  was  defined  based  upon  the  RSMs 

(Nos.  2,  4,  5  and  17) 

6  4  A  revised  estimate  ATj  =  120  zb  16  K  and  Uj  =  7.1  zb  1  m/s  obtained  using  the  RSMs 

7  2  A  revised  Tj—Too  =  115  K  was  selected  for  creation  of  the  RSMs  recognizing  that  the 

value  originally  used  (Tj  —  Too  =  66  K)  was  far  below  the  value  predicted  by  the  RSMs 

8  4,5  The  new  RSMs  yield  the  estimate  Tj  —  Too  =  105  zb  13  K  and  Uj  =  7.1  zb  1  m/s 

•  Result 

Quantity  Experiment  Predicted 

Tj  -  Too  107  zb  10  K  105  zb  13  K 
Uj  8.0  m/s  7.1  zb  1  m/s 
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Conclusions 


•  Developed  DDDAS  methodology  for  evaluation  of  fluid  thermal  systems 

—  Examples  are  optical  fibre  furnace  and  turbofan  combustor 

—  Need  for  complete  flowfield  simulation  to  optimize  system  performance 

—  Boundary  conditions  for  flowfield  simulation  are  not  completely  known  a  priori 

—  Non-intrusive  optical  measurements  (e.g.,  laser  diode  absorbance)  feasible  in  limited  region 

—  DDDAS  method  to  determine  complete  boundary  conditions  by  synergizing  experiment  and 
simulation 

•  Developed  DDDAS  method  to  determining  Tj  and  Uj 

•  DDDAS  method  predicts  Tj  —  Too  and  Uj  within  experimental  uncertainty 
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